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ABSTRACT 

We numerically investigate the possibility that a close stellar encounter could account for the high 
inclinations of the Kuiper belt, as originally proposed by Ida, Larwood and Burkert, however we consider 
encounters with pericenters within those explored by previous works. A star of mass Q.2Mq entering the 
solar system on a nearly parabolic, low inclination, retrograde orbit with perihelion of 50 AU can scatter 
~ 30% of the Kuiper belt into moderate inclination {i > 10°) and eccentricity orbits, while leaving the 
rest of the belt at low eccentricities and inclinations. This scenario provides a possible explanation for the 
dichotomy of the Kuiper belt, accounting for the scattered and classical populations, including objects 
with high eccentricities and inclinations. High eccentricity objects that were not placed into resonance 
with Neptune are more likely to be removed from the system by subsequent long timescale evolution. 
Consequently this scenario suggests that some Plutinos could have originally been objects in circular 
orbits which were scattered into resonant high eccentricity and high inclination orbits by the stellar 
encounter. Objects could have been placed into stable resonant regions with Neptune by the encounter. 
However long lived regions of dynamical stability should differ from those selected by resonance capture. 
This scenario predicts the presence of objects in resonances with Neptune that are unlikely to have been 
captured via Neptune's migration. 

Subject headings: Kuiper Belt — planetary systems: protoplanetary disks 



1. INTRODUCTION 

The Edgeworth Kuiper Belt (KB) is often discussed in 
terms of three dynamical families; 1) the classical Kuiper 
Belt objects which occupy semi-major axes 41 < a < 46 
AU and which have low eccentricities e < 0.25, 2) the 
resonant Kuiper Belt occupying the 3:2 and 2:1 mean mo- 
tion resonances with Neptune, denoted the Plutinos and 
Twotinos, respectively, and 3) the scattered Kuiper belt 
occupying extreme orbits with high eccentricities (e up to 
0.6) (e.g, Trujillo & Brown 2002). Most of the mass in the 
Kuiper belt appears to reside in the classical Kuiper belt 
(Luu & Jewitt 2002). All three dynamical families con- 
tain objects with surprisingly high inclinations (i > 15°; 
Trujillo & Brown 2002; Luu & Jewitt 2002; Brown 2001). 
Current compilations suggest a bimodal distribution, with 
~ 20% of the Kuiper belt objects in low-inclination or- 
bits of a few degrees and the remainder in highly inclined 
orbits with i > 20° (Brown 2001). Scenarios for the for- 
mation and evolution of the Kuiper Belt should explain 
the different Kuiper Belt Object (KBO) populations. 

Several mechanisms have been proposed to account for 
the high KBO inclinations, including scattering by Nep- 
tune during Neptune's migration (Gomes 2003; Levison 
& Morbidelli 2003), passage through vertical resonances 
during Neptune's migration or during the depletion of the 
solar nebula (Nagasawa & Ida 2000; Malhotra 1998), wave 
excitation (Ward 2003), perturbations from passing stars 



(Ida et al. 2000; Kobayashi & Ida 2001), and scattering 
by planets and planetesimals (Thommes et al. 2002; Petit 
et al. 1999; Kenyon 2002; Brunini & Melita 2002). While 
long term evolution can remove KBOs with low inclina- 
tions, it cannot explain the existence of classical KBOs 
with extremely high inclinations of 30° (Kuchner et al. 
2002). Current surveys appear to rule out the possibility 
of a distant planetary embryo (Morbidelli et al. 2002). 

In this paper we focus on the possibility that perturba- 
tions from a stellar encounter could account for the high 
KBO inclinations, as proposed by Ida et al. (2000). The 
probability of a stellar encounter that passes with a few 
hundred AU of the Sun at its current position is extremely 
unlikely. However, most stars are formed in stellar groups 
or clusters (Lada & Lada 2003), associated for 10 — 10^ 
years until they are tidally disrupted by molecular clouds 
and structure in the Milky Way disk (Lada & Lada 2003; 
Terlevich 1987; Theuns 1992; Bergond et al. 2001). Young 
stellar groups and clusters, can have stellar densities up 
to 10^ times higher than the Solar neighborhood (Lada 
& Lada 2003; Hillendbrand & Hartmann 1998). The Sun 
could have experienced a close stellar encounter when it 
was still in its birth cluster, (Ida et al. 2000; Adams & 
Laughlin 2001; Gaidos 1995). Other disk systems such 
as Beta Pictorus and HD 100546 probably have suffered 
close stellar encounters (Larwood & Kalas 2001; Quillen & 
Thorndike 2002). Alternatively, the Sun could have been 
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part of a wide binary which was later disrupted (e.g., as 
Furlan et al. 2003). 

As shown by Larwood & Kalas (2001) in their study 
of the Beta Pictorus system, and others (e.g. Kobayashi 
& Ida 2001) a close stellar encounter can scatter plan- 
etesimals. For an equal mass parabolic encounter with 
pericenter distance q, planetesimals with semimajor axis 
a > q/2 are scattered, whereas those with a < q/2 are uni- 
formly perturbed into higher eccentricity and inclination 
orbits (Larwood & Kalas 2001; Kobayashi & Ida 2001). 
This division depends on the timescale of the encounter 
compared to the orbital rotation period of the planetesi- 
mals. In the inner disk, the planetesimals rotate quickly 
compared to the encounter time and the collision can be 
treated with a secular approximation, averaging over the 
entire orbit, as done in the appendix by Kobayashi & Ida 
(2001). However, objects with larger semimajor axes are 
moving at speeds similar to the encountering star (assum- 
ing a parabolic orbit). Consequently the resulting per- 
turbation on outer planetesimals is strongly dependent on 
their longitudes, as shown explicitly by Larwood & Kalas 
(2001). If the collision is very fast compared to the ro- 
tational timescale of the planetesimals, as would be the 
case if the star was on a highly hyperbolic orbit, then the 
impulse approximation can be used to estimate the pertur- 
bations caused by the encounter (e.g., Eggers & Woolfson 
1996; Pfalzer 2003). 

Previous work considering the effect of a stellar en- 
counter on our solar system has considered parabolic or- 
bits with pericenter of size 100-200 AU (Ida et al. 2000; 
Kobayashi & Ida 2001). At this distance the velocity of 
the perturber at pericenter is similar to the rotational ve- 
locity of KBOs. For such an encounter the Kuiper belt 
is in a transition region. The KBOs experience increases 
in mean eccentricity and inclination which are dependent 
on their initial semimajor axes. Because the final KBO 
orbits in this region also depend on the initial longitude, 
the eccentricity and inclination distribution broadens fol- 
lowing the encounter. Ida et al. (2000); Kobayashi & Ida 
(2001) restricted their exploration to fairly large stellar 
pericenters q > 100 AU, which minimized the eccentricity 
variations in the KBOs caused by the encounter. Ida et 
al. (2000); Kobayashi & Ida (2001) suggested that Nep- 
tune's migration followed the encounter. KBOs with ini- 
tial eccentricities above 0.15 are unlikely to be captured in 
the 3:2 mean motion resonance with a migrating Neptune 
(Malhotra 1995; Hahn & Malhotra 1999). The encounters 
considered by Ida et al. (2000); Kobayashi & Ida (2001) 
were consistent with a scenario that included migration of 
Neptune since the encounters did not induce high KBO 
eccentricities. 

In this paper we extend the previous explorations by Ida 
et al. (2000); Kobayashi & Ida (2001) to include stars on 
orbits with smaller pericenters. In this regime the plan- 
etesimal perturbations are strongly dependent on the plan- 
etesimal longitude during the encounter. Consequently 
we expect a larger spread in the resulting planetesimal 
planetesimal orbital properties, as seen in the outer parts 
of the planetesimal disks simulated by Larwood & Kalas 
(2001); Ida ct al. (2000); Kobayashi & Ida (2001). Here 
we explore a regime where the Kuiper Belt is more effec- 
tively scattered by the encounter than previously works 



have considered. 

2. CONSTRAINTS ON PARAMETER SPACE OF THE 
ENCOUNTER 

In this section we discuss preliminary constraints on the 
orbit, mass and initial velocity of the stellar impactor. We 
maintain our hypothesis that the Kuiper Belt is strongly 
scattered by the encounter, but that it must retain low ec- 
centricity and inclination objects as are found in the clas- 
sical Kuiper Belt. Furthermore we assert that Uranus's 
orbit is not strongly perturbed by the encounter. 

The Sun/stellar encounter is described by 5 free param- 
eters, the mass of the star compared to the Sun, M^,/Mq, 
the velocity of the star distant from the Sun, Voo , the dis- 
tance of closest approach from the Sun or perihelion of the 
star, (/, the inclination of the star's orbit with respect to 
the ecliptic, i*, and the argument of perihelion, w,. For 
a pictorial view of these angles see Fig. 1 which is based 
on Fig. 1 by Kobayashi & Ida (2001). The velocity of 
the incoming star is related to the eccentricity of the orbit 
e*, where = (e* - 1) ^2iM@±MA^ and e* > 1 for a 

hyperbolic orbit. 

For a star to scatter a planetesimal effectively, the star's 
velocity must be greater than the circular velocity of the 
planetesimal. The velocity of a star in a hyperbolic orbit 
at perihelion is 

.g = (l + e.) ^^^^+^-) . (1) 

We compare the star's velocity at perihelion with the ve- 
locity, Vc, of a particle in a circular orbit with semi-major 
axis a; Vc = ^ GMqI a. Allowing Neptune to be in the 
scattered region, we require Vc ^ Vq at Neptune's a = 30 
AU which implies that 



ql^JTTT, < 30AU. (2) 

This sets pericenter q < 50AU if the star is on a parabolic 
orbit and if Neptune is in the scattered region. 

We now constrain e* based on the most probably stel- 
lar incoming velocity. Collisions are unlikely from field 
stars, so we expect an incoming stellar velocity of order a 
few km/s, typical of the velocity dispersion in young star 
clusters which have not yet dispersed or disrupted. 

/G(M0+M* 



1) 



(3) 



V 9 

The circular velocity is 4.2 km/s at a semi-major axis of 
50 AU, of order that which might be found in a young 
unbound cluster. Since we don't expect to be greater 
than a few km/s if the perturber is in a birth cluster, we 
infer that the encounter was likely to be nearly parabolic, 
(e* < 2). A hyperbolic orbit with a small pericenter from a 
faster field star would be possible, but it is much less likely 
than a parabolic one from a star in a dense birth cluster. 
If the encounter was a result of an interaction with a bi- 
nary, then the velocity of the encounter would depend on 
the nature of the binary and how it was disrupted. 

We can set a condition on the orbit of the impactor by 
specifying that Uranus does not significantly increase in in- 
clination or eccentricity following the encounter. Uranus 
is at small enough semimajor axis that we may assume 
that it moves quickly compared to the timescale of the en- 
counter. To estimate the affect of the encounter on Uranus 
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we use equations 15 and 16 from Kobayashi & Ida (2001) 
which are appropriate for parabolic encounters in a secular 
approximation. 

3/2 

«~ — 1= , ' I - I |sin2u| (4) 



8^2 v'l + M*/M0 \q, 

Based on Uranus's incUnation, 0.77°, and semimajor axis, 
19.2 AU, and assuming an encounter with q 50AU we 
find 

sin 2z* < 0.087. (5) 

If the star was low mass M* ~ 0.1, then a reasonably 
large range of inclinations is allowed. Since Uranus's ec- 
centricity increase is oc (a/q)^^^ (Kobayashi & Ida 2001), 
Uranus's low eccentricity is less of a constraint than its 
low inclination. If Uranus was not strongly perturbed by 
the encounter, then the stellar impactor was likely to be a 
low mass star or on a low inclination orbit. 

For particles with semi-major axes similar to the peri- 
ccntcr distance of the encounter, the perturbation caused 
by the encounter is strongly dependent upon the particle's 
longitude. Since these particles move slower than the star, 
we can use the impulse approximation to estimate the size 
of the perturbations. Each particle undergoes a change in 
its velocity of order 

bVg 

where b is tiic distance between the; object and the star's 
closest approach to it, and Vq is the; velocity of the 
star at pericenter (e.g., Binney & Tremaine 1987; Eg- 
gers & Woolfson 1996). The velocity change is in the 
direction toward the star at its position of closest ap- 
proach. Using cylindrical coordinates {R, z) where R is 
in the ecliptic, the position of the star's pericenter is 
R = q\/ cos^ ti)^ + sin^ tj* cos^ , z = gsinw* sini*. For 
an object opposite pericenter (180° away) 
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To allow the inclination of particles opposite pericenter to 
remain low, we restrict the encounter to low arguments of 
perihelion and low inclinations. 

3. SIMULATIONS OF HYPERBOLIC ENCOUNTERS 

Prom the previous section we found that if Neptune is in 
the scattered region, then q < 50AU. If the incoming star 
is part of the Sun's birth cluster then we expect e* < 2. To 
maintain Uranus's low eccentricity and inclination, we re- 
quire low inclination orbits and low mass stars. To ensure 
that the Kuiper Belt retains a low eccentricity and incli- 
nation population we require orbits with low inclinations 
and low arguments of perihelion. We take these estimates 
as guidelines for our numerical integrations. 

To investigate the actual affect of a such a close stellar 
encounter on a primordial cold Kuiper belt, we numer- 
ically integrated the two body Sun/stellar gravitational 
system. Our code is a conventional Burlisch-Stoer numer- 
ical scheme which considers only the force from Newtonian 
gravity. Simultaneously during the encounter we integrate 
1000 low mass particles (planetesimals) per run which were 

^ http://cfa-www.harvard.edu/iau/lists/TNOs.html 



initially placed in circular orbits in the ecliptic. These 
particles are influenced by the gravitational force from the 
Sun and star but do not themselves act on the Sun or star. 
The integrations were begun with the star at 500 AU from 
the Sun, and ended when the star reached an equivalent 
distance from the Sun following the encounter. We began 
our simulations with planetesimals distributed randomly 
in semi-major axis between 30 and 70 AU and randomly 
chosen mean longitudes. Following the encounter, the fi- 
nal semi-major axes, eccentricities and inclinations of the 
particles were tabulated. 

Based on our constraints from the previous section, we 
explored planetesimal eccentricity and inclination distri- 
butions for encounters spanning 0.1 < M^/Mq < 0.5, 
40 < g < 70 AU with differing orbital inclinations and 
arguments of perihelion (i*,a;*). Since we expect that the 
encounter could have been in the Sun's birth cluster, we 
restrict the incoming stellar velocity v^c to less than a few 
km/s, equivalent to nearly parabolic encounters, e* < 2. 
In this region of parameter space, we searched for final 
planetesimal distributions which contained both scattered 
and classical type KBO population analogs. 

In Fig. 2 we show the planetesimal distribution follow- 
ing a retrograde encounter with a star of M^/Mq = 0.2, 
e* = 1.4, pericenter q = 40AU, inchnation i» = 170° and 
argument of perihelion = 20°. In Fig. 2, histograms 
are also shown for all particles with final semi-major axes 
between 33 and 50 AU. Examination of this figure shows 
that the inclination and eccentricity distribution follow- 
ing the encounter is wider at larger semi-major axes as 
expected from a comparison between the timescale of the 
encounter with the orbital period of the objects. We see 
that a significant fraction, but not all, of the belt is scat- 
tered, indicating that the final eccentricity and inclination 
of an object depends on its longitude during the collision 
(as shown explicitly by Larwood & Kalas 2001). About 
30% of objects with final semimajor axes between 33 and 
50 AU have inclinations above 10° and eccentricities above 
0.1. However, low eccentricity (e < 0.05) and low inclina- 
tion (i < 6°) objects remain in the planetesimal distribu- 
tion. 

Fig. 2 can be compared to the distribution of known 
Trans-Neptunian objects (as listed by the Minor Planet 
Center Jan. 7, 2004 ^) which is given in the same for- 
mat in Fig. 3. To ensure accurate orbits, only objects 
that have been observed over more than one opposition 
have been plotted. We must keep in mind that the known 
population of KBOs is likely to be biased by the systemat- 
ics of the observations (e.g., as discussed by Bernstein et 
al. 2004; Brown 2001). The KBO has also had ~ 4 billion 
years to evolve dynamically following the hypothetical col- 
lision considered here (e.g., as explored by Kuchner et al. 
2002). 

While the simulation of the retrograde encounter shown 
in Fig. 2 does not exactly match the known KBO pop- 
ulation (the most obvious mismatch being the radial dis- 
tribution), it successfully predicts the combined presence 
of objects with high inclinations and eccentricities while 
retaining a population with low inclinations and eccen- 
tricities, as required by the known KBO distribution. The 
inclination and eccentricity distribution following the en- 
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counter have high inchnation and eecentricity tails, similar 
to those seen in the KBO comparison histograms. An ob- 
ject at the location of Neptune (30 AU) can remain in a low 
eccentricity and inclination orbit, so Neptune itself could 
have survived the encounter with minimal changes in its 
orbital elements. 

The distribution of orbital elements for the retrograde 
encounter can be contrasted, for example, with the pro- 
grade encounter shown in Fig. 4 which does not exhibit 
as broad a distribution of eccentricities and inclinations. 
The simulation (Fig. 2) we chose to directly compare to 
the current KBO population (Fig. 3) corresponded to 
a low inclination retrograde encounter. A prograde en- 
counter can more effectively perturb the disk because the 
collision effectively lasts longer in the frame rotating with 
each planetesimal. We also find prograde collisions which 
are capable of scattering the disk, and an example of one 
is shown in Fig. 4. While low inclination retrograde en- 
counters fail to produce very low inclination objects in the 
belt, we find that the low inclination prograde close en- 
counters fail to produce low eccentricity objects, requiring 
subsequent evolution to account for them, including the 
low eccentricity of Neptune itself. In addition, we have 
had difficulty identifying a region of parameter space in 
which a low eccentricity and inclination population is re- 
tained but a a significant fraction of the belt is scattered. 

We consider what effects subsequent evolution could 
have on the resulting planetesimal distribution shown in 
Fig. 2 for the retrograde encounter. Scattered objects at 
high eccentricities that were not placed in resonances with 
Neptune would be quickly removed from the system be- 
cause they cross Neptune's orbit and so could be scattered 
by Neptune. High ccc;entric;ity objcx'ts at the location of 
Pluto (40AU) would be removed, leaving a population of 
Plutinos which include high eccentricity and inclination 
objects. Alternatively because there remains a population 
of low eccentricity objects, subsequent migration by Nep- 
tune can resonantly capture objects into the 3:2 resonance 
accounting for the Plutinos. Because the low eccentricity 
objects include objects at high inclinations, the resulting 
captured Plutinos would also include high inclination ob- 
jects. 

The simulation shown in Fig. 2 does not necessarily 
require Neptune to migrate to account for the Plutinos 
since non resonant high eccentricity objects at 40 AU will 
be removed by subsequent scattering by Neptune. If we 
consider models which include Neptune's migration, then 
a closer (smaller pericenter) encounter than simulated in 
Fig. 2 would be required. This follows because low eccen- 
tricity objects captured at smaller semi-major axes must 
include high inclination objects such as are found past 40 
AU in this simulation. 

While no objects are found at extremely low inclinations 
in the retrograde simulation, (Fig. 2), subsequent evolu- 
tion may broaden the inclination distribution, reducing 
the inclinations of some objects. We began the simula- 
tions with zero inclination and eccentricity objects. Had 
we started the simulation with a distribution of inclina- 
tions and eccentricities, the final distribution would have 
included lower inclination objects. If the stellar orbit has a 
higher inclination than the simulation shown in 2, then fol- 
lowing the encounter the mean inclination is higher, how- 



ever more objects remain at low eccentricities, and the 
eccentricity distribution is narrower. 

Deep observational searches suggest that the Kuiper 
Belt is truncated at 50 AU (Allen et al. 2001) and does not 
contain a large population of "cold" or near zero eccentric- 
ity objects (Allen et al. 2002). Our simulated encounters 
do not truncate the disk, but instead increase the radial 
distribution of any previously defined edge (e.g. as shown 
by (Larwood & Kalas 2001) in similar simulations for the 
Beta Pic disk). It is possible that an edge in the KBO dis- 
tribution existed prior to the stellar encounter. To inves- 
tigate the effect of this edge, we ran encounter simulations 
with a a predefined edge at 45 AU in the planetesimal dis- 
tribution prior to the encounter. Following the encounter, 
the planetesimal distribution is shown in Fig. 5. The 
simulated encounter was retrograde with = O.25M0, 
= 1.3, u = 170°, = 20°, and go = 40AU. The 
resulting planetesimal distribution exhibits a rise in both 
mean eccentricity and inclinations with semi-major axis 
for objects past 45 AU. A rise in mean eccentricity is ex- 
hibited by the known KBOs shown in Fig. 3, however the 
known KBO distribution lacks an increase in mean incli- 
nation in the region same region. Consequently the ob- 
served increase in mean eccentricity with semimajor axis 
seen in the KBO population is unlikely to be explained by 
an encounter perturbing a pre-existing disk edge. Since 
the brightness of an object depends on its distance from 
the Sun to the fourth power, high eccentricity KBO's arc 
more likely to be discovered when they are near perihelion, 
even taking into account the time they spend at small dis- 
tances from the Sun compared to that spent at large dis- 
tances. This selection effect may have reduced the number 
of known low eccentricity KBOs objects at large semima- 
jor axis. 

4. DISCUSSION AND SUMMARY 

In this paper we have numerically investigated the ef- 
fect of a close stellar encounter on the primitive Kuiper 
Belt. We have extended our study to collisions with lower 
mass stars and smaller pericentcrs than have the pioneer- 
ing studies of Ida et al. (2000); Kobayashi & Ida (2001). 
We find that low mass, low inclination encounters can 
produce a dichotomy in the resulting planetesimal distri- 
bution, containing both scattered objects with high ec- 
centricities or high inclinations as well as low inclination 
and low eccentricity objects, similar to what is seen in the 
Kuiper Belt. In particular a O.2M0 star passing with a 
pericenter of 45 AU from the Sun, in a retrograde, highly 
inclined orbit produces inclination and eccentricity distri- 
butions with high end tails, similar to those seen in the 
KBO distributions. 

To more thoroughly test such a scenario it is necessary 
to numerically evolve the scattered population for the life- 
time of the solar system as done for the classical KBO 
population by Kuchner et al. (2002). High eccentricity 
objects at 40AU which were not placed in resonance with 
Neptune are likely to be removed from the Kuiper Belt 
by scattering with Neptune. Consequently the 3:2 reso- 
nance represents a region of stability; objects placed in 
this resonance could stay there for the lifetime of the solar 
system. Objects placed in this resonance by the encounter 
could explain the high inclination Plutino population. Ob- 
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jects currently residing in resonances that have low cap- 
ture probabilities such as the 5:2 resonance (Chiang et al. 
2003) might be more naturally explained with a scattering 
scenario rather than a Neptune migration scenario. Al- 
ternately one could regard the existence of such objects 
as a consequence or prediction of the impact scenario ex- 
plored here. Since long lived regions of dynamical stability 
differ from those selected by resonance capture, study of 
future measurements of the of KBO orbital distribution 
may support or reject a scenario which includes a close 
stellar encounter. 

Low eccentricity objects remain in the belt following 
the encounter, which implies that there remains a pop- 
ulation which can be captured into resonance by subse- 
quent migration by Neptune (e.g., Nagasawa & Ida 2000), 
however Neptune's migration may not be needed to ex- 
plain the highest Plutino eccentricities. This may reduce 
the extent of migration required by recent migration mod- 
els (Gomes 2003; Levison & Morbidclh 2003). Hahn & 
Malhotra (2003) have suggested that the resonant KBO 
populations might be better explained by Neptune's mi- 
gration into a previously "hot" Kuiper Belt, a situation 
which could occur after the encounter we have simulated 
here. 

Following the stellar encounter, we expect many objects 
to be removed from the belt by dynamical evolution, sug- 
gesting that there was a period of time when short period 
comets were produced in great numbers. One interesting 
possibility is that this epoch corresponded to the period 
of "Late Heavy Bombardment" when the solar system was 
a few xlO^ years old. This possibility would require a 
cometary origin for the lunar impacts associated with the 
Late Heavy Bombardment Era (for further discussion on 
the cometary vs. asteroidal origin see Levison et al. 2000, 
2001). Because of the large mass of material likely to be 
sent into the inner solar system by Neptune, this epoch is 
also very likely to be the time during which Neptune mi- 
grated. Because of the narrowness of Neptune's mean mo- 
tion resonances, a close stellar encounter probably did not 
occur after the migration of Neptune. This follows because 
such an encounter would have effectively removed objects 
from resonance. The encounter explored here would also 
place objects in the 2:1 resonance with Neptune, and so 
may not naturally account for the low fraction of Twotinos 
compared to Plutinos unless there was a pre-existing edge 
to the Kuiper Belt. 

While a close encounter with a low mass star could ac- 
count for the inclination distribution and dichotomy of or- 
bital properties seen in the Kuiper Belt, as seen in our 
simulations, it does not truncate the disk. Because they 
pull tidal arms outward, stellar encounters do not tend to 
sharply truncate disks (e.g. Pfalzer 2003). If the Kuiper 
Belt contains an edge at 50 AU as suggested by (Allen et 
al. 2001), then it is much more likely to have been caused 
by a companion, either stellar or planetary. We note that 
the encounters simulated here did remove ~ 10% of the 
planetesimals in the region 30-70AU objects from the So- 
lar system. The encounter could have ejected a primordial 
planet at a semimajor axis of ^ 50 AU, leaving only its 
signature behind in the KBO distribution. 

If a stellar encounter did perturb the young Kuiper 
Belt, then the limits on the likely stellar birth cluster of 



the Sun can be somewhat relaxed from the lower limits 
placed by Adams & Laughlin (2001), as suggested by (de 
la Fuente Marcos & de la Fuente Marcos 2001). In partic- 
ular a higher stellar density birth cluster, similar to that 
in the Orion nebula with a few times 10^ stars pc~^ (Hil- 
lendbrand & Hartmann 1998), or a longer lasting open 
cluster, would be required to make the hypothetical col- 
lision explored here probable during the lifetime of the 
cluster. The probability P, of an object of cross section 
a undergoing a collision while moving through a stellar 
population of number density n, with velocity dispersion 
V ^ 1 km/s during time t is P = nvat (e.g., Adams & 
Laughlin 2001). Due to gravitational focusing, the impact 
parameter b, of an encounter is larger than the pericen- 
ter distance q. For an encounter with a small pericenter 



(when G(Mq + M» ) / q > Voo), the impact parameter is 

6 ^ VG(MQ+M,)/g 
9 ■ 



For an encounter with the Sun with 

pericenter g ~ 50 AU, and Voo = Ikm/s, we estimate a 
cross section a = 7r6^ ^ 3 x 10"^ AU^. In the Smi's birth 
cluster the probability that such an encounter occurred 
during time t 



0.16 
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t 



,3 X lO^AUV VllOMyr 
(8) 

We therefore find that a collision proposed here would 

have been only a moderately improbable event (15 per- 
cent level) if the Sun were born in a dense cluster similar 
to that in the Orion nebula. Here we have only considered 
the probability of a close collision and have not restricted 
the orientation of the encounter with respect to the eclip- 
tic. Had we taken into account the narrow region of al- 
lowed orbital inclinations and arguments of perihelion, our 
proposed encounter would be about 100 times less prob- 
able. The formation time of the KBOs (> 10® years) is 
longer than that used in the above estimate, (Kenyon & 
Luu 1998; Stern & ColweU 1997), suggesting that the Sun 
resided in a longer lived and more massive birth cluster, 
such as an open cluster. A longer timescale would al- 
low the close encounter to occur after the formation of 
the KBOs. Alternatively the Sun could have been part 
of a wide binary, and the close encounter could have oc- 
curred when the binary was disrupted (e.g., one of the 
components of T Tauri has just recently been disrupted, 
see Furlan et al. 2003). 

Currently there are number of examples of young sys- 
tems which show clear evidence of external perturbers, 
including the AeBe stars HD 141569A which has spiral 
structure and asymmetries excited by its binary compan- 
ion HD 141569B,C (Augereau & Papaloizou 2003; Quillen 
et al. 2004), and HD 100546 which exhibits spiral struc- 
ture which could have been excited by star at a pericenter 
of about 500 AU (QuiUen et al. 2004; Grady et al. 2001). 
The tilt in Bet Pictorus's disk has been explained with an 
stellar encounter with a similar sized peric;enter (Larwood 
& Kalas 2001). Structure seen in submillimeter emission 
of older systems have recently been interpreted in terms of 
high eccentricity outer planets. For example models for the 
submillimeter distribution for Vega include an outer high 
eccentricity Jupiter mass planet (Wilner et al. 2002) and 
that for Epsilon Eridani include a Neptune mass moder- 
ate eccentricity planet (Quillen & Thorndike 2002). While 
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the hypothetical encounter with our solar system explored 
here did not significantly increase Neptune's eccentricity, 
it easily could have, had the incoming star approached the 
Sun from a different orientation. This suggests a possible 
scenario which could account for the observed morphology 
of the the dusty debris disks which require outer eccentric 
planets. Future study may determine where and when 
close stellar encounters are likely to affect young solar sys- 
tems as well as test and further explore the possibility that 
such an encounter occurred in our own Solar system. 

We thank S. Thorndike, E. Chiang, E. Ford, P. Varniere, 
and A. Frank for helpful discussions. This material is 
based on upon work supported by the National Aeronau- 
tics and Space Administration under Grant No.xxx is- 
sued through the Origins of Solar Systems program. This 
research was supported in part by the National Science 
Foundation to the KITP under Grant No. PHY99-07949. 
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Fig. 1.— 



Orbital angles for parabolic and hyperbolic encounters, based on the similaj: figure by Kobayashi & Ida (2001). 
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Fig. 2. — Planetesimal distribution following a retrograde hyperbolic encounter with a star. The star has M* = 0.2Mq, orbit pericenter at 
q = 40AU, inclination i* = 170°, argument of perhelion = 20° and eccentricity e« = 1.4. The lower left shows the eccentricity distribution 
as a function of semi-major axis. The upper left shows the inclination distribution as a function of semi-major axis. The right hand side 
shows histograms for planetesimals with semi-major axes between 35 and 50 AU. This low inclination encounter scatters a significant fraction 
of the planestesimals, while retaining a low inclination and low eccentricity analog to the classical Kuiper Belt population. 




Fig. 3. — Known Trans-Ncptunian objects from the minor planet center database as of Jan 7, 2004. To ensure accurate orbits, only 
objects that have been observed over more than one opposition have been plotted. The known KBO distribution is displayed identically to 
the simulation shown in Figure 2. The known KBO distribution of orbital elements has high eccentricity and inclination tails similar to those 
seen in the simulation presented in Fig. 2. The known KBO distribution is affected by observational selection affects and has had 4 billion 
years of dynamical evolution since formation. High eccentricity objects that are not in mean motion resonances with Neptune are likely to 
have been removed from the Kuiper belt. 
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a (AU) N(33<a<50) 

Fk;. 4. Similar to Figure 2 except the encounter is prograde. This simulation corresponds to an encounter with Af, = O.IMq, e, = 1.5, 
it = 20°, LOt = 20°, and go = 60AU. Prograde encounters can also cause a dichotomy in the scattered population similar to what is seen in 
the Kuiper Belt. However, we have failed to identify a region in parameter space that retains low eccentricity objects as well as low inclination 
ones, and still scatters significant fraction of the belt. 




Fig. 5. — Similar to Figure 2 except the planetesimal disk prior to the encounter was truncated at 45 AU prior to the encounter. This 
simulation corresponds to an encounter with M* = O.25M0, e* = 1.3, i* = 170°, u!» = 20°, and qo = 40AU. Following the encounter, objects 
initially near the edge of the disk can be scattered to higher eccentricity, and semi-major axis orbits. Past 45 AU, both the mean eccentricity 
and inclination increase with semi-major axis. A previously truncated disk cannot explain the eccentricity increase with semimajor axis seen 
in the Kuiper belt because there is no increase in mean inclination past 45 AU in the Kuiper belt. 



